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The physics during the inflationary stage of the universe is of quantum nature involving extremely 
high energy densities. Moreover, it is out of equilibrium on a fastly expanding dynamical geometry. 
We complement here the 1999 Chalong e Lectures on out of equilibrium fields in self-consistent 
inflationary dynamics [ ^stro-ph / 0006446 investigating inflation driven by the evolution of highly 
excited quantum states. These states are characterized by a non-perturbatively large number of 
quanta in a band of momenta and with zero or nonzero expectation value of the inflaton scalar 
field. They represent the situation in which initially a non-perturbatively large energy density is 
localized in a band of high energy quantum modes and are coined tsunami- waves. The self-consistent 
evolution of this quantum state and the scale factor is studied analytically and numerically. It is 
shown that the time evolution of these quantum states lead to two consecutive stages of inflation 
under conditions that are the quantum analogue of slow-roll. The evolution of the scale factor during 
the first stage has new features that are characteristic of the quantum state. During this initial stage 
the quantum fiuctuations in the highly excited band build up an effective homogeneous condensate 
with a non-perturbatively large amplitude as a consequence of the large number of quanta. The 
second stage of inflation is similar to the usual classical chaotic scenario but driven by this effective 
condensate. The excited quantum modes are already superhorizon in the first stage and do not 
affect the power spectrum of scalar perturbations. Thus, this tsunami quantum state provides a 
field theoretical justification for chaotic scenarios driven by a classical homogeneous scalar field of 
large amplitude. 
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I. INTRODUCTION 

A wealth of observational evidence from the temperature anisotropics in the cosmic microwave background 
strongly points towards inflation as the mechanism to produce the primordial density perturbations |^]. 
Thus, inflationary cosmology emerges as the basic theoretical framework to explain not only the long- 
standing shortcomings of standard big bang cosmology but also to provide a testable paradigm for structure 
formation^-|Q. The recent explosion in the quantity and quality of data on temperature anisotropics 
elevates inflation to the realm of an experimentally testable scenario that leads to robust predictions that 
withstand detailed scrutiny ^. 

However at the level of implementation of an inflationary proposal, the situation is much less satisfactory. 
There are very many different models for inflation motivated by particle physics and most if not all of them 
invoke one or several scalar flelds, the inflaton(s), whose dynamical evolution in a scalar potential leads to 
an inflationary epoch[^-0. The inflaton fleld is a scalar fleld that provides an effective description for the 
flelds in the grand uniflcd theories. Furthermore there is the tantalizing prospect of learning some aspects 
of the inflationary potential (at least the part of the potential associated with the last few e-folds) through 
the temperature anisotropics of the cosmic microwave background]^. 

Most treatments of inflation study the evolution of the inflaton via the classical equations of motion for 
the chosen scalar potential instead of making a quantum field treatment of the dynamics. That is, the 
effect of quantum fluctuations is neglected in the dynamics of the inflaton. Furthermore, since inflation 
redshifts inhomogeneities very fast, the classical evolution is studied in terms of a homogeneous classical 
scalar field. The quantum fleld theory interpretation is that this classical, homogeneous fleld configuration 
is the expectation value of a quantum field operator in a translational invariant quantum state. While the 
evolution of this coherent field configuration (the expectation value or order parameter) is studied via classical 
equations of motion, quantum fluctuations of the scalar field around this expectation value are treated in 
a linear approximation for the high wavenumbers that lead to the seeds for scalar density perturbations of 
the metric [^-[0. The large amplitude modes that dominate the energy of the universe during inflation are 
mimic by the classical homogeneous field. Moreover, the slow roll approximation is used in all cases. 

A fairly broad catalog of infiationary models based on scalar field dynamics labels these either as 'small 
field' or 'large field' [^. In the 'small field' category the scalar field begins its evolution with an initial 
value very near the origin of the scalar potential and rolls down towards larger values, an example is new 
inflation|g|, U, In the 'large field' category, the scalar field begins very high up in the potential hill and 
rolls down towards smaller values, an example is chaotic inflation ^ ||. 

It is only recently that the quantum dynamics of the scalar flelds in the coupled evolution of matter 
and geometry has been studied self-consistently[^, |T^. This work is associated with the dynamics of non- 
equilibrium phase transitions in models that fall, broadly, in the 'small fleld' category. 

This subject is reviewed in the 1999 Chalonge Lectures on out of equilibrium flelds in self-consistent 
inflationary dynamics. The conclusion of these works is that a treatment of the quantum fluctuations that 
couple self-consistently to the dynamics of the metric provides a solid quantum fleld theoretical framework 
that justifles microscopically the picture based on classical inflation. At the same time these studies provide a 
deeper understanding of the quantum as well as classical aspects of inflation and inflationary perturbations. 
They clearly reveal the classicalization of initial quantum fluctuations!^ [To| , and furnish a microscopic 
explanation (and derivation) of the effective, homogeneous classical inflaton ||10|| . 

Very recently, a quantum dynamical treatment of models whose classical counterpart are large field models 
is proposed in ref . . The classical description in these models begins with a homogeneous infiaton scalar 
with very large amplitude 4> ^^p;[§, @' 01' i-^' "^^^V ^^S^ up in the scalar potential well. 

The initial quantum state may be mixed (described by a density matrix in Fock space) or a pure state 
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(described by a vector in Fock space). The initial state is characterized by the expectation value of the 
inflaton (order parameter) and the spectrum of excitations (initial particle distribution). The order parameter 
corresponds to the classical inflaton field in the classical limit whereas the initial particle distribution describes 
the spectrum of excitations in the initial state. 

It is not obvious whether an initial particle distribution can give rise to inflation. As shownin ref. 
there exists classes of initial particle distributions leading to efficient inflation. In particular, initial states 
with zero order parameter and a band of excited modes can lead to an inflationary epoch. This extends the 
set of possible initial states that leads to chaotic inflation making it a more natural description of the early 
universe. 

In the customary treatment of chaotic inflation (classical chaotic inflation 0) all the energy is contained 
in the classical (space independent) field, so the field modes do not contribute neither to the energy nor to the 
background dynamics. This is clearly a very special choice of initial conditions. We shall call this scenario 
classical chaotic inflation. 

We review in these lectures the dynamics that results from the evolution of a quantum state which drives 
the dynamics of the scale factor through the expectation value of the energy momentum tensor The 
main idea behind this approach is akin to the experimental situation in ultrarelativistic heavy ion collisions, 
wherein an initial highly excited state (heavy ions with very large energy) lead to the formation of a plasma 
that expands and cools |^|. Recently this situation has been modeled by considering the evolution of highly 
excited quantum states coined tsunami- waves in[|l3|, p^ . 

The goals of the Tsunami Inflation: 

The ideas and concepts in refs.||l^, |^ |l^ were adapted in ref.||ll| to study the self-consistent dynamics 
of the metric and the evolution of a highly excited quantum state with the goal of providing a quantum 
description of large field inflationary models without assuming an expectation value for the scalar field. Novel 
quantum states were considered, which are the cosmological counterpart of the tsunami-waves introduced 
inO, |l^, |l^ with the following properties [pT|: 

ya. quantum state with non-perturbatively large number of quanta in a momentum band, i.e, a large 
number of high energy excitations. We first solve explicitly the case of a narrow band and then we consider 
more general distributions of quanta. 

a) vanishing expectation value of the scalar field. 

The rationale behind considering these quantum states is that they provide a natural description of a 
situation in which a state of large energy density composed mainly of energetic particles evolves in time. 
The main idea is to model the situation conjectured to drive inflationary cosmology: that is an initial state 
with a large energy density and pressure which in turn couples to the metric leading to inflation. Efficient 
inflation should follow independently of the details of the initial state. In our case, this means that efficient 
inflation is to be obtained for different shapes of the initial particle distribution. That is, we want to avoid 
a fine tuning of the initial conditions. 

[In the usual studies (classical chaotic inflation) all of the initial energy is in the classical zero mode while 
the quantum fluctuations are taken to be in the ground state]. 

After inflation, energy transfer from the heavier to lighter particles results in reheating and particle 
production that eventually excites the light sector and leads to a radiation or matter dominated phase. 

The conditions under which a Tsunami quantum state leads to inflationary dynamics is established and 
the self-consistent evolution of this quantum state and the space-time metric is studied in detailpT|. 

We emphasize that we are not proposing here yet a new model of inflation. Instead we focus on inflation 
driven by the evolution of a quantum state, within the framework of familiar models based on scalar flelds 
with typical quartic potentials. This is in contrast with the usual approach in which the dynamics is driven 
by the evolution of a homogeneous classical field of large amplitude. 

Brief summary: We find that inflation occurs under fairly general conditions that are the quantum equiv- 
alent of slow-roll. There are two consecutive but distinct inflationary stages: the first one is completely 
determined by the quantum features of the state. Even when the expectation value of the scalar field van- 
ishes at all times in this quantum state, the dynamics of the first stage gives rise to the emergence of an 
effective classical homogeneous condensate. The amplitude of the effective condensate is non-perturbatively 
large as a consequence of the non-perturbatively large number of quanta in the band of excited wavevectors. 
The second stage is similar to the familiar classical chaotic scenario, and can be interpreted as being driven 
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by the dynamics of the effective homogeneous condensate If . The band of excited quantum modes, if not 



superhorizon initially they cross the horizon during the first stage of inflation, hence they do not modify 
the power spectrum of scalar density perturbations on wavelengths that are of cosmological relevance today. 
Actually, in the explicit examples worked out in ref.[|ll], the excited modes are initially superhorizon due 
to the generalized slow-roll condition. Therefore, in a very well defined manner, tsunami quantum states 
provide a quantum field theoretical justification, a microscopic basis, for chaotic inflation, explaining the 
classical dynamics of the homogeneous scalar field. 

In section II we introduce the quantum state, obtain the renormalized equations of motion for the self- 
consistent evolution of the quantum state and the scale factor. In section III we provide detailed analytic and 
numerical studies of the evolution and highlight the different inflationary stages. In section IV we discuss 
generalized scenarios. The summary of results is presented in the conclusions. An appendix is devoted to 
technical details and the equations of motion for mixed states. 



II. INITIAL STATE AND EQUATIONS OF MOTION 

As emphasized in the introduction, while most works on inflation treat the dynamics of the inflaton field 
at the classical level, we use a quantum description of the inflaton. 

We focus on the possibility of inflation through the dynamical quantum evolution of a highly excited initial 
state with large energy density. Consistently with inflation at a scale well below the Planck energy, we treat 
the inflaton field describing the matter as a quantum field whereas gravity is treated semiclassically. 

The dynamics of the classical space-time metric is determined by the Einstein equations with a source 
term given by the expectation value of the energy momentum tensor of the quantum inflaton field. The 
quantum field evolution is calculated in the resulting metric. 

Hence we solve self- consistently the coupled evolution equations for the classical metric and the quantum 
infiaton field. 

We assume that the universe is homogeneous, isotropic and spatially flat, thus it is described by the metric, 

ds^ = dt^ - a^{t) dx^ . (2.1) 

Anticipating the need for a non-perturbative treatment of the evolution of the quantum state, we consider 
an inflaton model with an A'^-component scalar inflaton field with quartic self-coupling. We then invoke 
the large N limit as a non-perturbative tool to study the dynamics H, |l^ |6j. This choice is not only 

motivated by the necessity of a consistent non-perturbative treatment but also because any grand unified 
field theory will contain a large number of scalar fields, thus justifying a large N limit on physical grounds. 

The matter action and Lagrangian density are given by 

= / d^x£„, = [ d-^xa^it) 



(x) - - 
^ ' 2 



(V$(x-))' 



- ^('i>(^)) 



(2.2) 



vm = ^ $2 + A ($2^ ' -I- i e 7^ $\ (2.3) 

and we will consider > 0, postponing the discussion of the case < to future work. Here Tl{t) stands 
for the scalar curvature 



The ^-coupling of $^(a;) to the scalar curvature Tl{t) has been included in the Lagrangian since it is necessary 
for the renormalizability of the theory. 
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The discussion of the alternative inflationary mechanism that we are proposing and the physical description 
of the quantum states becomes more clear in conformal time 



T = 



dt' 



in terms of which the metric is conformal to that in Minkowski space-time 

ds^ = a^{T) {dT^ - rfx^) . 
We introduce the conformally rescaled field 

*(T,x)=a(i)%x) 
in terms of which the matter action becomes 

S[^] = J dTd'x niidr^r - (V*)^] - a\T)V ^ 



(2.5) 

(2.6) 
(2.7) 

(2.8) 



Since we are interested in describing the time evolution of an initial quantum state, we pass on to the 
Hamiltonian description in the Schrodingcr representation. This procedure begins by obtaining the canonical 
momentum conjugate to the quantum field, n(T, x), and the Hamiltonian density H{T,x) 



n(T,x) = *'(T,x), 



H(r,x) = in2 + i(v*)2 + a4(r)y 



H{T) = J d^yiHiT, 



a(T) 



a^(r)|*2 



(2.9) 



where the prime denotes derivative with respect to the conformal time T. 
In the Schrodinger representation the canonical momentum is given by 



n''(r,x) = -i 



(5*''(T,x) 



; a=l,...,N. 



(2.10) 



The time evolution of the wave-functional T 

. d 



H 



is obtained from the functional Schrodinger equation 

d 



(2.11) 



The implementation of the large N limit begins by writing the field as follows 

*(x, T) = (a(x, T), 7r(x, T)) 

= (x/]VV(T)+x(x,T),fr(x,T)), 



(2.12) 



where we choose the 'I'-axis in the direction of the expectation value of the field and we collectively denote 
by TT the N — 1 perpendicular directions. 

V(r) = (a(x, T)) 
(7r(x, T)) = (x(x, T)) = , (2.13) 

where the expectations value above are obtained in the state represented by the wave- functional T 
introduced above. 
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The leading order in the large N limit can be efficiently obtained by functional methods (sec refs.[^, |T^, 
|r^ , |l5| , |l6|| and references therein). The contributions of x to the equations of motion are subleading (of 
order 1/iV) in the large N limit [|, |l6). 

It is convenient to introduce the spatial Fourier modes of the quantum field 



(fx 7f(x, T) 



(2.14) 



In leading order in the large TV limit, the explicit form of the Hamiltonian is given 

2 



Id, 14 la If 



H{T) = NVKi{T)-^ (y^{A-^^^)\ +Y,H,{T) 

\ k / k 



Hk{T) ^ 



1 



2 Jtt/j • (57r_; 



a\T) 



M\T) 



n{T) 
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A V' A (7r2) 



2 a2(r) 2 iVa2(r) 



(2.15) 
(2.16) 
(2.17) 



where V is the comoving volume. We assume spherically sy mme tric distributions in momentum space. 

That is, in the large N limit the Hamiltonian operator ( ^.9[ ) becomes a time dependent c-number con- 
tribution plus a quantum mechanical contribution, Hk{T), given by a collection of harmonic oscillators 
with time -dep en dent frequencies, coupled only through the quantum fluctuations (tt/c • 7r_fe). 

In eqs.( 2.15 )-( ^.17 ) the scale factor a(T) is determined self-consistently by the Einstein-Friedmann equa- 
tions. 



A. Tsunami initial states 



To highlight the main aspects of the inflationary scenario proposed here, and to establish a clear difference 
with the conventional models, we first focus our discussion on the case of vanishing expectation value, i.e, 
'0(T) = 0, and a pure quantum state. The most general cases with mixed states described by density matrices 
and non-vanishing expectation value of the field are discussed in detail in sections IIIC and IV and in the 
appendix. 

For ip{T) = the quantum Hamiltonian ( ^.15 ) becomes a sum of harmonic oscillators with time dependent 
frequencies that depend on the quantum fiuctuations. Therefore, we propose a Gaussian wave- functional of 
the form 



T 



:AAx(r)n 



(2.18) 



The functional Schrodinger equation (2.11) in this case leads to evolution equations for the normalization 
factor AfriX) ^''^^ the covariance kernel Ak(T) whose general form is found in the appendix (see also Q). 
The evolution of the normalization factor is determined by that of Ak, while the equation for is 



^A',iT)^A 



(2.19) 



where primes refer to derivatives with respect to conformal time. As described in the appendix for the 
general case, the above equation can be linearized by defining (see appendix) 



(2.20) 
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where the mode functions ipk satisfy the equation 

In terms of these mode functions the self-consistent expectation vahie 



(2.21) 



N 



1 

N 



N 



(Pk 
(2^ 

N 



Wk\ 



(2.22) 



We now must provide initial conditions on the wave functional to completely specify the dynamics. Choosing 
the initial (conformal) time at T = with a(T = 0) = 1, the initial state is completely specified by furnishing 
the real and imaginary parts of the covariance at the initial time. We parameterize these as |17| 



AR^k{0)^^k ; Ai^k{0)^u;k{0)Sk 
Choosing the Wronskian of the mode functions (fikC^) and its complex conjugate to be 

'Pk ^k -^'kf*k = 2i 

determines the following initial condition on the mode functions (see appendix) 

1 



fkiO) = 



^'k{0) ^ - K(0)4 + i ^k] ^k{0) ■ 



(2.23) 



(2.24) 



(2.25) 



where we also used eqs.(2.2C) and (2.23). 

An important alternative interpretation of these mode functions is that they form a basis for expanding 
the Heisenberg field operators (solution of the Heisenberg equations of motion) 



tt{x, T) = 



(2^ 



(2.26) 



with ak]a\, annihilati on an d creation operators, respectively, with canonical commutation relations. The 
Wronskian condition ( 2.24 ) ensures that the 7r(x, T) fields and their conjugate momenta obey the canonical 
commutation relations at equal conformal times. 

The physical interpretation of these initial states is highlighted by focusing on the occupation number of 
adiabatic states as well as on the probability distribution of field configurations. 

• Occupation number: it is at this point where the description in terms of conformal time proves to be 
valuable. In conformal time the Hamiltonian in the large N limit is that of a collection of harmonic 
oscillators with time dependent frequencies. It is then convenient to introduce the adiabatic occupation 
number operator 



nk{T) = 1 



L^k{T) 



(2.27) 



with Hk given by eq. (2.15). 

In particular the occupation number at the initial time is given by (see appendix) 



nk = {nk{0)) = 



K(0)-»fc]Vc.^(0) SI 

4 wfc(o) nk 



(2.28) 



Here, the special case with Vlk = ujkiO) and 6k = corresponds to the adiabatic vacuum (rik = ). 
Instead, we study an initial state in which a band of wave- vectors are populated with a non-perturbatively 
large number of particles. More precisely, we consider initial states for which 



rik = 0[j 



inside the excited band, 
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rik = outside the excited band. 



(2.29) 



where A is the quartic self-couphng. 
This is accomphshed by choosing, 



1 
1 



O 



1 



Acjfc(O) 



1 



and 



inside the excited band, 
Sk = outside the excited band. 



(2.30) 



These initial states are highly excited, the expectation value of the energy momentum tensor in these 
states leads to an energy density ^ 1/A and are, therefore, non-perturbative. We will refer to the case 
where the excited band is narrow as the narrow tsunami. 

It must be stressed that the particle distribution Uk alone partially determines the initial state. As we 
see from eq.( p.25 ), the initial state is completely defined specifying two functions of fc : Qk and Sk ■ 



• Probability distribution: an alternative interpretation of these initial states is obtained by focusing on 
the probability distribution of field configurations at the initial time. It is given by 



V[7f] 



T 



vE-iT = 



:AA(o)n 



-ill, •?fc-7f_fc 



(2.31) 



An intuitive quantum mechanical picture of the wave-functional for the modes in the excited band 
is the following. At the initial time the instantaneous Hamiltonian corresponds to a set of harmonic 
oscillators of frequencies uJk{Q), while the width in field space of the initial Gaussian state is determined 

by ri^^^^. For a mode in the vacuum state l/y/Tlk ~ 1/ ^/lujJO) and the typical amplitudes of the field 
are tt^ ~ 1/ y^ujk{0) which is the typical width of the potential well. While for a mode inside the excited 
band the width in field space is l/\/TTk ^ l/y^X ujk{0) [see eq.( 2.30 )]. Thus, large amplitude field 
configurations with nk^ko ^ l/-\/A ujk{0) ^ l/^/uJk{0) have a probabihty of 0(1), i.e, large amplitude 
configurations within the band of excited wave-vectors are not suppressed. That is, the width of the 
probability distribution for these modes is much larger than the typical size of the potential well and 
there is a non-negligible probability for finding field configurations with large amplitudes of 0(1/A). 

These highly excited initial states had been previously proposed as models to describe the initial stages 
of ultrarelativistic heavy ion collisions and had been coined 'tsunami waves '[[I |l|, H). They represent 
spherical shells (in momentum space) with large occupation numbers of quanta, describing a state with 
a large energy density with particles of a given momentum. 

To highlight the difference between these states and other proposals for infiation, we now summarize the 
important aspects of these tsunami-wave states before we study the conditions under which such states lead 
to inflation. 



Properties of Tsunami states 



Pure states: the states under consideration, defined by the wave-functional of the form ( 2.18 ) are pure 
states unlike some other proposals that invoke a thermal distribution of particles, hence a mixed state 
density matrix as initial state. 

Vanishing expectation value of the scalar field: Unlike most models of classical chaotic inflation in 
which the scalar field obtains an expectation value, taken as a classical fiel d, th e expectation value 
of the scalar field in the tsunami- wave states given by the wave-functional ( [2.18 ) vanishes. We also 



consider tsunami states with a non-zero expectation value for the scalar field [see sec. IV]. 
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• Highly excited initia l mode s: the tsunami-wave state described by eq. (2.18) with the covariance kernel 
given by eqs. ( ^.23 ), ( 2.30 ) describes a state for which the modes inside the excited band are occupied 
with a non-perturbatively large 0{1/X) number of (adiabatic) quanta. We remark that very high energy 
modes, those that will become superhorizon during the last 10 or so e-folds, hence are of cosmological 
importance today, must be in the vacuum state so as not to lead to a large amplitude of scalar density 
perturbations 1^, 0, ||. 

This type of quantum states is clearly a novel concept, it presents an alternative to typical inflationary 
scenarios that invoke the dynamics of a classical scalar field which in most cases ignore or bypass the 
quantum dynamics. 

The tsunami-wave states described above are the simplest states and will be the focus of our study. These 
states can be generalized to describe mixed-state density matrices and to also allow for an expectation value 
of the scalar field. These generalizations are described in sections III.C and IV and in the appendix and are 
found to lead qualitatively the same features revealed by the simpler pure states. 



B. Back to comoving time: renormalized equations of motion 



Having set up the initial value problem in terms of the tsunami-wave initial wave-functionals, the dynamics 
is now complet ely d etermined by the set of mode equations eqs. (2. 21) with eqs.(2.16)-( p.l7 ) and the initial 
conditions eqs. (2.25). It is convenient to re- w rite the equations of motion in comoving time. This is accom- 
plished by the field rescaling given by eq.(2.7) which at the level of mode function results in introducing the 
comoving time mode functions fk (t) related to the conformal time ones tpk (T) as 



fk{t) 



Mr) 

a{t) 



The equations of motion in comoving time for these mode functions are 



h{t)+iH{t)h{t) 



a^{t) 



-M\t) 



fk{t) = 



A 



d^k 



\Ut)? 



fk{0) 



/fc(0) 



(27r)3 

[uJk{0)dk + HiO) + ink] /fc(0) . 



The Einstein-Friedmann equation are. 



o(t) 



3M2, 



Po = (Too) , 



(2.32) 

(2.33) 
(2.34) 
(2.35) 

(2.36) 



where the expectation value is taken in the time evolved quantum state. It is straightforward to see that 
the expectation value of the energy-momentum tensor has the perfect fluid form, as a consequence of the 
homogeneity and isotropy of the s y stem P, p^ . 

Thus the set of equations (2.33)-(2.36) provide a closed set of self-consistent equation for the dynamics of 
the quantum state and the space-time metric. 



Renormalized Equations of Motion in the Large N limit 



The set of equations that determine the dynamics of the quantum state and the scale factor need to be 
renormalized. The fleld quantum fluctuations 



N 



d^k 
2(27r)= 



\fk{t)\' , 
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requires subtractions which are absorbed in a renormalization of the mass, coupling to the Ricci scalar and 
coupling constant. The expectation value of the stress tensor also requires subtractions (but not multi- 
plicative renormalization). Since the divergence structure is determined by the large energy, short distance 
behavior, the band of excited modes does not influence the renormalization aspects. Therefore, we use the 
extensive work on the renormalization program which is available in the literature referring the reader to 
references^, ^ for details. We here summarize the aspects that are most relevant for the present discussion. 
First, it is convenient to introduce the following dimensionless quantities, 

H{t) k 
HV = ; g = — ; 



m t 



, — ; ^q = — ; 9 = ^ 

m m 57r^ 

where m and A stand for the renormalized mass of the inflaton and the renormalized self-coupling, 
respectively In terms of these dimensionless quantities we now introduce the dimensionless and fully 
renormalized expectation value of the self-consistent field as 



fqij) = Vm fk{t) 



(2.37) 



S(r) 



2m2 
q^dq 



1 



q a(T)2 



e(g-l) /Af(r) 

2q^ \ m? Qm? 



(2.38) 



where the terms subtracted inside the integrand renormalize the mass, coupling to the Ricci scalar and 
the coupling constant |lOj. The dimensionless and renormalized expressions for the energy density e and 
pressure p are given by 



e(r)^ 



A 



2N m4 



)R = 



2 

p{t) = 



+ 1 / q^dq{\f,{T)\^-S,{q,T) + 



a2(T) 



[|/,(r)p-^2(g,r)] 



A 



2N m4 



{p + e)iT)=g / q'dq\\f,{T)\^~S,{q,T) + 



3a2(r) 



[|/,(r)p-52(g,r)] 



(2.39) 



Where the renormalization subtractions 5*1 and 5*2 are given by,|^, ^] 



Siiq, t) 



S2iq,T) 



1 



+ 



a4(T) 2(7a4(r) 

e{q-r 



[5(7 



8g3 a^(T) 
1 



1 



a{TfB{T) + 3a{T)d{T)B{T) - 4d^{T)B{T) 



qa?'{T) 2q^a?{T) 
B{t) = a2(T)[l + gS(r)] . 



B{t) 



e(g-i) 

8q^ a^{T) 



3B(t)^ +a(T)— a(T)B(T 
dr L 



(2.40) 



We choose here ^ = (minimal coupling), the renormalization point k = |r7i| and a(0) = 1. 

In summary, the set of coupled, self-consistent equations of motion for the quantum state and the scale 
factor are 

^2 



1 + .gS(r) 



Mr) = 




a2(r) 

/,(0)--K S, + h{0) + tn,]f,{0) 



5^(0) 



no) 

6to2 



(2.41) 
(2.42) 

(2.43) 
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plus the Einstein-Friedmann equation of motion for the scale factor 



h\T)=L^e(T) 



where 



16 TriV m? 



(2.44) 



with 5S(t) and e(T) given by eqs. (2.38) and (2.39) respectively. 

In order to implement the numerical analysis of the set of eqs. (2.41)-(2.42), (2.38) and ( ^.44 ) we introduce 
an ultraviolet momentum cutoff A. For the cases considered here we choose A ~ 200 and found almost no 
dependence on the cutoff for larger values. As befits a scalar inflationary model, the scalar self-coupling 
is constrained by the amplitude of scalar density perturbations to be A '--^ 10"^^ 0, ^, ^ implying that 
g < 10"^'^. Therefore, the subtractions can be neglected because Si '-^ 0(gA**) < 10 . 

The initial state is defined by specifying the fig and Sq . We determine the range of these parameters Jl, 
and Sq by the excitation spectrum for the tsunami-wave initial state, as well as the condition that lead to 
infiationary stage. This will be studied in detail in the next section. 



III. TSUNAMI INFLATION 



As emphasized in the previous section, the scenario under consideration is very different from the popular 
treatments of inflation based on the evolution of classical scalar infiaton field|^, |[ ||, ^, In these scenarios 
all of the initial energy is assumed to be in a zero mode (or order parameter) at the beginning of inflation and 
the quantum fluctuations are taken to be perturbatively small with a negligible contribution to the energy 
density and the evolution of the scale factor. 

In contrast to this description, our proposal highlights the dynamics of the quantum states as the driving 
mechanism for inflation. The initial quantum states under consideration correspond to a band of quantum 
modes in excited states, thus the name 'tsunami- wave' ]T^, [T^ , |5). This initial state models a cosmological 
initial condition in which the energy density is non-perturbatively large, but concentrated in the quanta 
rather than in a zero mode. 

We now study under which general conditions such a state can lead to a period of inflation that satisfles 
the cosmological constraints for solving the horizon and entropy problems entailing the necessity for about 
60 e-folds of inflation. 

It is understood that inflation takes place whenever the expansion of the universe accelerates, i.e, 

a 

-^h^ + h = - — [e + 3p]>0, (3.1) 



with L given in eq. (2.44) and e and p given by eqs. (2.39) 



While our full analysis rely on the numerical integration of the above set of equations, much we learn by 



considering the narrow tsunami case 



A. Analytical study: the narrow tsunami case 



Before proceeding to a full numerical study of the equations of motion, we want to obtain an analytic 
estimate of the conditions under which a tsunami initial quantum state would lead to inflation. 

Our main criterion for such initial state to represent high energy excitations is that the number of quanta 
in the band of excited modes is of 0(1/A). This criterion, as explained above, is tantamount to requiring that 
field configurations with non-perturbative amplitudes have non-negligible functional probability. Progress 
can be made analytically by focusing on the case in which the band of excited field modes is narrow i.e, 
its width A/c is such that A/c <C or in terms of dimensionless quantities Aq/qo <^ 1. We introduce the 
following smooth distribution 



(3.2) 
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with LOq given by eq. ( 2.43| ) and Afn a normalization constant that fixes the value of the total energy. 
In addition, we choose 5q — ~h{Q)/ujq as we discuss below in eq.(3.13). 

This initial distribution posses the main features of the tsunami state described in the previous section. 
Since g ^ 1, we have for q ^ qo, 



1 1 
> 1 



a 



9 



i»l 
9 



corresponding to highly excited states. While for — go| ^ 



1 



1 



(3.3) 



(3.4) 



Thus, these modes are in a quantum state near the conformal (adiabatic) vacuum at the initial time, with 
Uq the number of quanta defined by eq. ( p. 2^ ) in terms of dimensionless variables. For these distrib ution s 
(narrow tsunamis), the integral over mode functions for the quantum fluctuations 5S(r) [given by eq. (2.38)] 
is dominated by the narrow band of excited states with mode amplitudes ~ 1/ ^ and can be approximated 

by 



5E(t) - g Aq ql |/,„(t)|2 + 0[g) + 0[g Aq) ^ |0,„(t)|2 
where we have introduced the effective q^ mode 

(t^qoij) = \/9 go /goM- 



(3.5) 



(3.6) 



we note that the initial condition ( 2.42| ) and the tsunami-wave condition (3.3) entail that despite the presence 
of the coupling constant in its definition, the amplitude of the effective go niode is of 0(1). 
The equation of motion for the effective go-mode takes the form 



4>qa (t) + 3 /i(t) (t) -f 

The scale factor follows from 



h\T) = e{T) , 



+ l + l0«oWI 



(3.7) 



(3.8) 



with energy and pressure, 
e(r) 
(p + e)(r) 



l0,„(T)r 



4 



(3.9) 



where we have neglected terms of 0{g). 

We will refer to the set of evolution equations ( |3.7D -(3.£) as the one mode approximation evolution equa- 
tions. 

In particular, within this one-mode approximation, the acceleration of the scale factor obeys 



a(r) 
a{r) 



41 



Therefore, the condition for an inflationary epoch, a > 0, becomes 



(3.10) 



(3.11) 
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A sufficient criterion that guarantees inflation is the tsunami slow roll condition 



(3.12) 



The initial conditions ( 2.42| ) and the condition that Qq^ ^ g ^ 1 imply that the tsunami slow roll condition 
( 3.12 ) at early times is guaranteed if 6qg is such that 



\uJqJq„+h{0)\^l (3.13) 

Hence tsunami- wave initial states that satisfy the tsunami slow- roll condition ( 3.13| ) lead to an inflationary 
stage. 



Moreover, in order to have slow roll ( 3.12 ) at later times, the effe ctive friction coefficient 3/i(r) should be 
larger than the square of the frequency in the evolution equations (2.41). That is, 



2 



1 + .gS(0) 



3h{0) 



< 1 . 



(3.14) 



{i.e. the go-mode should b e dee p inside the overdamped oscillatory regime). Eq.(3.14) implies that h{0) 3> 1 
and this together with eq.( 3.13|) implies that Sqg must be negative. 

A remarkable aspect of the narrow tsunami state is that it leads to a dynamical evolution of the metric 
similar to that obtained in classical chaotic in flation ary scenarios in the slow roll approximation j^, ^. In 
particular the expression for the acceleration (3.10) and the tsunami slow roll condition (3.12) are indeed 
similar to those obtained in classical chaotic inflationary models driven by a homogeneous classical field (zero 
mode). However, despite the striking similarity with classical chaotic models , we haste to add that both the 
conditions that define a tsunami state and the tsunami slow roll condition ( 3.12| ) guaranteed by the initial 
value ( |3.13 ) is of purely quantum mechanical origin in contrast with the classical chaotic-slow- roll scenario. 
Furthermore, we recall that the expectation value of the scalar field vanishes in this state. 



Early time dynamics 



Under the assumption of a tsunami wave initial state and the tsunami slow- roll condition (3.12) the 
contribution (/)g„ (r) in the energy and in the pressure [see eq. ( |3.9| )] can be neglected provided, 



3a2(r) 



Mr)? 



(3.15) 



We call ta the time scale at which this rely no longer holds. Furthermore, we can approximate 4'qg (r) by 



TA < 



,(0) 



""90 



(0) 



(3.16) 



This condition is fulfilled at least for ta ^ 1 due to the tsunami slow-roll condition ( 3.12 ) . 
During this interval the Friedmann equation (3.8) takes the form. 



r«(r)i 


2 







2 1-^90(0)^ + ^10,0(0)1' + 



2a2(r) 



D 



E 



(3.17) 



where we used that gS(0) = |(/)5g(0)|2. This equation is valid as long as the characteristic time scale of 
variation of the metric is shorter than that of the mode (pq^ij). 
The preceding equation can be integrated with solution 



a(T) = sin\\(\fET 



a(r) 
a(r) 



E > 



14 



E coth 



h{r) 



E 



sinh^ (\fET + c 



(3.18) 



where the constants D, E and c are given by, 

E^L 



2 / 5^0 , 5^0 




(3.19) 



and gI]o = <?S(0) = |<^a„(0)p. 

We see from eqs. (3.18) that during this interval there is an inflationary stage with an accelerated expansion 

= E > 0. We also see that /i(r) decreases with time until it reaches the constant value \/E that 

rin 

deter mines the onset of a quasi-De Sitter inflationary stage. In ref. is shown that the solution in eqs. 
( 3.18 ) is valid for t < ta, where 



1 



ArgSinh 



VSgpE 





1 








ArgSinh ^ 




~7e 





'£V3^(l+.gSo/2) 
, V2(l+5So) 



(3.20) 



This initial inflationary period with a decreasing Hubble parameter exists provided the r.h.s. is here positive, 
i.e., 



V3^ 



(3.21) 



In order to distinguish this phase from the later stages, to be described below, we refer to this early time 
inflationary stage as 'tsunami-wave inflation' because the distinct evolution of the scale factor during this 
stage is consequence of the tsunami-wave properties. 
At r = we have: 



^o('^a) ^ (0) = Va^o 
/3^<7o 



a{TA) 



1 + .gSo 



Hta) ^ 



1 + gSo 

3h{TA) 



4>qo (ta) 



E 



(3.22) 



For T > TA, qo/a{T) << 1 and the physical wavevectors in the excited band have red-shifted so much that 
all terms containing go become negligible in the evolution equations. Therefore, all modes in the excited 
band evolve as an effective q = mode. Hence for t > ta the dynamics of the scale factor is described by an 
effective homogeneous zero mode and describes a different regime from the one studied above. Such regime 
is akin to the classical chaotic scenario. 



2. The effective classical chaotic inflationary epoch 

For T >> TA, when Qq/o^ ^ I "/"g P / 1 "/"g P ^ 1 tti6 physical momenta corresponding to the comoving 
wavevectors in the excited band have redshifted to become negligible in the equations of motion. The 
dynamics is now determined by the following set of equations for the effective zero mode and the scale 
factor, 

4>,jT)+3h{T) 0,„(r) + [l + |0,„(T)p]0,„(r) = 0, 

h^{T) ^ e{T) , (3.23) 
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where the energy and pressure are given by, 

1 , • 

ip + e){r) = |^,„(r)p 



(r)P + il0,o(T)P + 7l'^.o(r)|^ 



(3.24) 



The initial conditions on (j)q^-, and (l)qg are determined by their values at the time ta, while the slow-roll 
condition ( 3.12| ) determines that the imaginary parts of (f>qg and (pq^ are negligible. 

Therefore, after ta the dynamic is identical to that of a classical homogeneous field (zero mode) 



Veffir) = Re[(j)q„{T)] , 

that satisfies the equations of motion, 

iieff +3hf]eff + {l+Veff)Veff = 



(3.25) 



h'iT)=L^eiT) 



(3.26) 



with energy and pressure, 



e(r) 
(p + e)(T) 



1 2 ^1 4 



and initial conditions [using cq. ( 3.22 )] 



VeffiTA) 

Veffir a) 



Veff 



(pqoiTA) 



(3.27) 



'90 (0) = V 5^0 

1 + 5S0 



3h{TA) 



Veffir a) ■ 



(3.28) 



Where the value of VeffirA) is determined by the slow roll condition (^ (^^^(r^) ~ 0), the evolution eq. 
(3.7) and a(r^) and /i(ta) are given by eq. ( [3.22 ). 

When gTiQ <^ 1, the quadratic term in the potential dominates, and we can integrate the previous equations 
to obtain 



Veffir) = Veffir a) - |^('^ 



ta) , (for gSo < 1) • 



(3.29) 



This evolution is similar to that of classical chaotic inflationary models]^, ^, ^. Therefore for t > ta when 
the physical momenta in the excited band have redshifted so much that their contribution in the equations of 
motion of the quantum modes and the energy and pressure become negligible, the evolution of the quantum 
modes and the metric is akin to a classical chaotic inflationary scenario driven by a homogeneous c-number 
scalar field. This equivalence allows us to use the results obtained for classical chaotic inflation. Thus, as the 
classical slow roll condition (|77e//| ^ \Veff\) holds, the evolution of the effective scalar flcld is overdamped 
and the system enters a quasi-De Sitter inflationary epoch. This inflationary period ends when the slowly 
decreasing Hubble parameter becomes of the order of the inflaton mass, i.e, 3h ^ 1 + rjeff- ^^^^ stage the 
effective classical field exits the overdamped regime and starts to oscillate, the slow roll condition no longer 
holds and a matter dominated epoch (|?7e//| ~ l'7e//| =^ p ^ 0) follows. 



3. Number of e-folds 



An important cosmological quantity is the total number of e-folds during inflation. As discu s sed a bove, 
there are two different inflationary stages, the flrst one is determined by the equations (3.17)-(3.18) and 
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characterized by a rapid fall-off of the Hubble parameter approaching a quasi-De Sitter stage. This new 
stage has been referred to as the tsunami-wave inflationary stage above to emphasize that the dynamics is 
determined by the distinct characteristics of the tsunami-wave initial stage. 



The second stage is described by an effective zero mode and the evolution equations (3.26, 3.27) and is 
akin to the chaotic inflationary stage driven by a classical homogeneous scalar field. The crossover between 
the two regimes is determined by the time scale ta (in units of the inflaton mass) and given by eqs.(3.2C) 
at which the contribution from the term q^/a?[T) to the equations of motion becomes negligible. Therefore 
there are to distinct contributions to the total number of e- folds, which is given by 



N,{qo, h{0)) = \oga{TA) + iVe(0, /i(t^)) 



(3.30) 



where a(r/i) is given by eq. ( 3.28 ) and iVe(0, /i(t^)) is just the number of e-folds for classical chaotic inflation 
with an initial Hubble parameter hijA)- 

We can express h{TA) as a function of qo and h{0), 



Kta) 



The number of e-folds during the first stage, is given by 



log a(rA) ~ log 



I + 5S0 



(3.31) 



(3.32) 



The expression for the number of e-folds during the following, chaotic inflationary stage simplifies when 
(7E0 ^ 1. In this case the quadratic term in the potential dominates, and we can obtain simple analytical 
expressions 



iVe(0, h{TA)) - — Tlljf{TA) = — gSo - — 



2 l + qt 



1 1 



(for 5E0 < 1) 



(3.33) 



We see that the number of efolds grow when go decreases at fixed initial energy eg. That is, we have more 
efolds when the energy is concentrated at low momenta. 



4. In summary 



Before proceeding to a full numerical study of the evolution we summarize the main features of the 
dynamics gleaned from the narrow tsunami case to compare with the numerical results. 

• The conditions for tsunami- wave inflation are i) a band of excited states centered at a momentum ko 
with a non-perturbatively large 0{l/g) number of quanta in this band, and ii) the tsunami slow-roll 
condition eq .(|3.12| ). These conditions are guaranteed by the initial conditions on the m ode fu nctio ns 
given by eq.( |2.42| ) with the tsunami-wave distributions of the general form given by eqs. (3.2), (3.13). 



There are two successive inflationary periods. During the first one, described in sec. HIAl, the 
dynamics is completely characterized by the distinct features of the tsunami-wave initial state, the 
Hubble parameter falls off fast and reaches an approximately constant value ^/E that chara cterizes the 
quasi-De Sitter epoch of inflation of the second period. The second stage, described in sec. [II A 2 can 
be described in terms of an effective classical zero mode and the evolution of this effective mode and 
that of the Hubble parameter are akin to the standard chaotic inflationary scenario. 

The tsunami-wave initial state can be interpreted as a microscopic justiflcation of the classical chaotic 
scenario described by an effective classical zero mode of large amplitude. The amplitude of this effective 
zero mode is non-perturbative as a consequence of the non-perturbative 0{l/ X) number of quanta in the 
narrow band of excited modes. Thus the initial value of the effective, classical zero mode that describes 
the second, chaotic inffationary stage, is completely determined by the quantum initial state. 
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• An important point from the perspective of structure formation is that the band of excited wavevectors 
centered at qo either correspond to superhorizon modes initiaUy, or all of the excited modes cross the 
horizon during the first stage of inflation, i.e, during the tsunami stage. This is important because the 
chaotic second stage of inflation which dominates during a longer period guarantees that the band of 
excited modes have become superhorizon well before the last 10 e-folds of inflation and hence cannot 
affect the power spectrum of the temperature anisotropics in the CMB. The fact that the tsunami- wave 
initial state is such that the very high energy modes (necessarily trans-Planckian) that cross the horizon 
during the last 10 e-folds and are therefore of cosmological importance today are in their (conformal) 
vacuum state leads to the usual results from chaotic inflation for the power spectrum of scalar density 
perturbations. 

Although these conclusions are based on the narrow tsunami case, we will see below that a full numerical 
integration of the s elf-consistent set of equations of motion confirms this picture. 

In sections III C and |^ we show how this results can be easily extended to more general particle distri- 



butions and more general initial states. 



B. Numerical example 

To make contact with familiar models of inflation with an inflaton field with a mass near the grand 
unification scale, we choose the following values of the parameters: 



M 



= 10- 



A = 10 



-12 



iV = 20 



PI 



(3.34) 



where the number of scalar fields = 20 has been chosen as a generic representative of a grand unified 
quantum field theory. For these values we find 

^ = 3.35 . 10« . 

3 A 

As an example we shall consider an initial energy density po — (Too) — 10^^ Mpi. Thus, the initial value 
for the Hubble parameter is Ho = ^/S-Kpo/iMpi = 3.53 • 10^^ GeV{^ 1.654 • 10^^ km/s/Mpc). These initial 
conditions in dimensionless variables gi ve ep = 2.50 and h{Q) — 2890. 

In addition, the slow roll conditions ( |3.14 ) imply: 



3/i(0) 



< 1 



which in this case results in 



qa < 95 



We choose q^ — 80.0, and initial conditions in eq. (2.42) with D.q and 5q given by eq. (3.2) and (3.13). These 



initial conditions satisfy the tsunami slow roll condition. 



h{Q)\ < 1 



(3.35) 



Furthermore, we take Ag = .1 and Mn is ad juste d by fixing the value gS(0) ~ 5E0 which for the values 
chosen for eo and by eqs.(|2!39|) and (U) and ( |3.12[) gives gEo = 7.81 • 10"''. 



Figure |l| displays eg vs go alon g lines of cons tant number of e-folds , while figures display the solution 
of the full set of equations ( ^.4l| )-( p.44| ) with ( 2.3E ). An important feature that emerges from these figures 
is that for the set of parameters that are typical for infiationary scenarios and for large values of go = ko/rn 
(but well below the Planck scale) the number of e-folds obtained is more than sufficient as shown by fig.P. 
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FIG. 1: Tsunami inflation: isolines of constant number of efolds obtained from eq.(3.33) (valid for gT,o <^ 1), for 
m = IQ-^Mpi, X = 10"^^ and N = 20. 



We also show that the dynamics of the full set of equations ( 2.41 )-( 2.44 ) with ( 2.39 ) is correctly approx- 
i mat e d by the narrow tsunami case studied in the previous sub sectio n s: th e one mode approximation [eqs. 
(|3.5[)-(3.9)], the ear ly tim e analytical formulae (for r < ta) [eqs. ( 3.18 )-( 3.19 )], and the effective classical field 
(for T > Ta) [eqs. (13^2^)1. The agreement between the analytic treatment and the full numerical evolution 
is displayed in figures 

The early time analytic expressions predict an inflationary period during which the Hubble parameter falls 
off fairly fast, that lasts up to ta ^ 0.133 [eq. (3.20)] reaching an asymptotic value of h{TA) = 36.2 [eqs. 
(3.19) and (3.22)]. The one mode approximation gives the same prediction /i(ta) = 36.1, and numerically 
evolving the full set of equations we find /i(t^) = 35.5. Thus, we see from this values and from figs. || and ^ 
that both approximations are fairly accurate for early times. 

After TA, the geometry reaches a quasi-De Sitter epoch. We have shown in the previous subsection that 
after the time ta the evolution equations for the one mode approximation reduce to those of an effective 
classical field. The effective zero mode approximation correctly predicts the dynamics in this epoch as can 
be gleaned from figures 

While the stage of early tsunami infiation up to ta results in only 8.5 efolds, the following quasi-de 
Sitter stage described by the effective classical scalar field lasts for a total of 1900 efolds. For the values of 
param eters chosen above, gSg ^ 1, hence we can estimate the number of efolds with eq. (3.33). Using eq. 



( 3.30 ) we obtain a total of 1970 efolds while the one mode approximation yields 1960 efolds. Both results 
agree with the full numerical solution of the equations (see fig. 

Furthermore, as stated above infiation ends when h ^ — ^ 
follows. 



after which a matter dominated epoch 
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FIG. 2: Tsunami inflation: Early time /i(r). Ha = h{TA) is the asymptotic value for the early period that ends 
at TA ~ 0.133. For m = W'^Mpianck, A = 10"" and = 20. Initial conditions: po = 10"^M^(, go = 80.0 and 
Aq = O.lgo. 



Other distributions and other states: 



The validity of the physical picture that emerges from the previous a naly tic and numerical study is not 
restricted to pure states or narrow distributions of the form given by (3.2). We have also studied more 
general distributions and mixed states: 

Other distributions: 

The narrow tsunami case where a single quantum mode go dominates the dynamics has been extremely 
useful to study the dynamics in the previous section. The generalization to the case with continuous distri- 
butions of g-modes can be easily obtained making the changes: 

\Kir)\' 



9 J q'dq |/,(r)p , 
9 Jq'dq |/,(r)p , 
9 I q^dqq' |/,(r)|2 



(3.36) 



The two stages of inflation are always present for such continuous modes distribution as long as the following 
generalized slow-roll conditions is fulfilled 



9 q^dq |/,(r)p « g dq |/,(r)| 



(3.37) 
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FIG. 3: Tsunami inflation: /i(r) for r > ta- The early time analytic approximation gives Ha = 36.1 (also with the 
one mode approx.), numerically we obtain h{TA) = 35.5. Same parameters and initial conditions as in fig. S. 



that imposes on Sq the condition \u!q Sq + h(0) \ <C 1. 

The effective zero mode in the second stage of inflation is now given by 

^eff(^)^9 J dq\fq{T)\K 

Our numerical study with general distributions reveals that the analytical picture obtained by substituting 
eq.( 3.36| ) in sec. IIIA correctly reproduce the dynamics. 

Other (mixed) states: Although we have focused for simplicity on tsunami pure initial states, we have 
also investigated the possibility of mixed states. Mixed state density matrices and their time evolution are 
discussed in the appendix. The mixing can be parametri zed in terms of angles Qk as given in equation (A4) 
and the number of (conformal) quanta are given by eq. ( A18). The only relevant changes that occur are in 
the integrals for S(t) ; £(t) ; p{t) in which 



l/,(r)P-|/,(r)pcoth 



l/.(r)P-|/,(r)pcoth 



(3.38) 



Tsunami and slow-roll conditions on the mode functions given by eqs. ( 2.30 ) and ( |3.13 ) lead to tsunami- wave 
inflation followed by chaotic inflation just as discussed above. In the narrow tsunami case the only change 
is that the effective go-mode is rescaled by the mixing factor, i.e, 



.Mp->|0,oMl'coth 



l'^,oWP^|<^,oMPcoth 



e 



(3.39) 
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FIG. 4: Tsunami inflation: it shows that there is accelerated expansion (inflation) up to times r ~ 109. Same 

parameters and initial conditions as in flg. 

It is also illuminating to contrast the tsunami- wave mixed states with the more familiar thermal mixed states. 
The latter are obtained by the choice 

ng = u;g ; ^, = ; 6, = (3.40) 

with T some value of temperature. In this case it is straightforward to see that (quantum) equipartition 
results in that the contributions of the modes and their time derivatives to the energy and pressure are 
of the same order (|/<j(t)P ~ [^(0)^ + '^g] l/gl''')!^)- Hence, for these thermal mixed states the tsunami 
slow-roll condition is not fulfilled. This is obviously not surprising, such a choice of thermal mixed state 
leads to a FRW epoch which is not inflationary. Hence the tsunami-wave initial conditions along with the 
generalized slow-roll conditions lead to two successive inflationary epochs in striking contrast to the familiar 
mixed thermal states. 



IV. GENERALIZED CHAOTIC INFLATION 

The previous analysis, confirmed by the numerical evolution of the full self-consistent set of equations 
leads to one of the important conclusions that tsunami-wave initial states provide a microscopic justification 
of the chaotic inflationary scenario. 

We have focused our discussion on initial states with vanishing expectation value of the scalar field (order 
parameter) and where the energy is concentrated in a momentum band (tsunami initial states). This choice 
brings to the fore the striking contrast between this novel quantum state and the usual classical approach to 
chaotic inflation. 
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FIG. 5: Tsunami inflation: (/S(r), after ta ~ 0.133, it plays the role of an effective classical fleld. Same parameters 
and initial conditions as in fig. ^ 



A relevant question to pose is what would be the dynamics if the initial state aUows for a non-vanishing 
expectation value of the scalar field along with some of the initial energy localized in excited quanta. In this 
section we address this question and study in detail this alternative scenario, which we refer to as generalized 
chaotic inflation to distinguish from the tsunami- wave state studied above. This generalization thus includes 
both cases: the classical chaotic inflation in the limit when there are no excited modes, as well as the tsunami 
initial state when all of the energy is localized in a band of excited modes and the expectation value of the 
field vanishes. 

The relevant equations of motion in comoving time f or the mode functions in this case are given by ( |A2C| ) in 



the appendix. Along with the dimensionless variables (2.37) it is also convenient to introduce a dimensionless 
expectation value as 



A 

2m2 



(4.1) 



In this generalized case with rj the equ ations of motion for the mode functions /^(t) (in terms of 
dimensionless variables) are the same as in eq. ( 2.41 ) after the replacement ^^(t) g^{T) + ??^(t) and the 
equation of motion for i]{t) is given by 



V{0) = Vo 



3M^)^ + [l + ^'(r)+5S(r)]7?(7 

; ^7(0) = m 



(4.2) 
(4.3) 



The Einstein- Friedmann equation is given by ( p. 44 ) but with the energy density and pressure now given by 



e(r) 



(5S 
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FIG. 6: ln[a(r)] vs. r. Same parameters and initial conditions as in fig. 



(p + e)(r) = i,^+g I dqi^\f,\^-S^{q,T) + £^[\f,\^~S2iq,r)]^ 



(4.4) 
(4.5) 



where the renormahzation subtractions Si ; 5*2 are obtained from those given by eq. ( 2.40 ) upon the 
replacement gYl — > + 77^ . 

From this expression we see that for fixed (large) energy density Status: RO 

there are two different possibilities: if the zero mode squared 77^(0) is larger than the quantum fluctuations 
gT, and g J q'^ dq \fq\'^, the dynamics is basically similar to that in the usual chaotic inflationary scenarios. 
This corresponds to most of the initial energy density to be in the zero mode and little energy density in 
the band of excited states. On the other hand, for small t]^{Q) most of the initial energy density is in the 
tsunami quantum state and the initial dynamics is akin to the rj = case. To quantify this statement and 
clarify the interplay and crossover of behaviors between the rj = and the generalized chaotic case, we now 
resort again to the narrow tsunami case, which highlights the essential physics. 

The relevant equations are: i) the equations of motion for the effective qo~iTiode (|3. 



(j)q„{T) +3/i(t) (r) + 
and for the zero mode r/ 



4 



l + r,2(r) + |0,„(r)|^ 



(4.6) 



fi(T) + 3/i(t),7(t) + [1 + r;2(r) + |0,„(T)n rKr) = 



(4.7) 
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FIG. 7: Tsunami inflation: p{T)/e{T). It shows the onset of a matter dominated epoch after the quasi-De Sitter 
stage. Same parameters and initial conditions as in fig. S. 



and the Hubble parameter given by (3.S) with the energy density given by 



IKir^- (4.8) 



The acceleration of the scale factor is now given by 



a{r) _ 2 



(4.9) 



where again we have neglected t he re normalization contributions and terms of 0{g) consistently in the weak 
coupling limit g ^ 1. From eq. (4.9) the generalized condition for an inflationary epoch (within the narrow 
tsunami case) becomes 



which is fulfilled if the following generalized slow roll condition holds 



'-"qo 



{r)\^+v'iT)«\Kir)\'+v'ir) 



(4.10) 



(4.11) 



Under these conditions and from cqs.( |4.6| )-(4.7) and the dynamics of the scale factor driven by the energy 
density eq.( [l.^ ) we can now distinguish the following different inflationary scenarios. 

• Tsunami dominated: When (qq + 1) |0go(O)p ^ V^iO) the excited states in the tsunami-wave carry 
most of the initial energy density. In this case the results of the previous section apply and the scale 
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factor takes the form as in eq.( |3.17 ) with D given by eg. ( [3. 19 ), and E given by eq.(3.19) but with 



gTiQ gYiQ + 7y (0). T here are two consecutive inflationary stages as in the previous section. The first 



described by eq.( 3.18D , lasts up to the time scale ta defined by 



4 



ia^{TA) 



at which the red-shift of the momentum go is such that qf)/a{TA) <^ 1. The secondary stage is a usual 
classical chaotic inflationary epoch determined by the dynamics of an effective zero mode given by 



(4.12) 



because for r > ta, qQ/a?{T) <^ 1 and the effective equation of motion for (j>qg{T) is the same as that 
for ry(r). 



Zero mode dominated: When ?7^(0) 3> (f^o + 1) I'Pg 



the energy density stored in the zero mode is 



much larger tha n that contributed by the excited states in the tsunami-wave. In this case the energy 
density eq.(4.8) is completely dominated by the zero mode. The ensuing dynamics is the familiar 
classical chaotic scen ario driven by a classical zero mode, without an early stage in which the scale 
factor is given by eq.(3.18) which is the hallmark of the tsunami- wave dynamics. 



This analysis in the narrow tsunami case does highlight the important aspects of the dynamics in a clear 
manner, allowing a clean separation of the two cases described above. We have carried a full numerical 
integration of the equations of motion that reproduce the results described above. The criterion for the 
crossover between tsunami-wave and classical chaotic inflation is determined by the relative contributions 
to the energy density from the quantum fluctuations in the tsunami wave state as compared to the energy 
density of the zero mode . 



The previous results [eqs.(4.1)-(4.12)] can be easily generalized for generi c con tinuous distributions of 
modes and for mixed states. One has just to make the changes indicated in eq.(3.36) for generic distributions 
and in eq.( 3.38| ) for mixed states. 

The generalized slow-roll conditions takes then the form: 



7?2(r)+5 / q^dq\fgiT)\' coth 



«'7'(t)+.9 / q^dq\fgiT)\^ coth 



during the first stage of inflation. 

The dynamics is tsunami dominated provided. 



g / q' dq{l + q')\f,m" coth 



e„ 



> T,2 



The effective zero mode in the second stage of inflation is now given by 



Vtffir) ^V\r)+g / q' dq |/,(t)|^ coth 



These results have been verified by numerical integration of the full set of evolution equations (2.39)-(p^44) 



V. CONCLUSIONS 



We have presented here infiation in typical scalar field theories as a consequence of the time evolution of 
a novel quantum state. This quantum state is characterized by a vanishing expectation value of the scalar 
field, i.e, a vanishing zero mode, but a non-perturbatively large number of quanta in a momentum band, 
thus its name-tsunami-wave state. 
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This state leads to a non-perturbatively large energy density which is localized in the band of excited 
quantum modes. We find that the self-consistent equations for the evolution of this quantum state and the 
scale factor lead to inflation under conditions that are the quantum analog of slow-roll. 

The self-consistent evolution was studied analytically and numerically in a wide range of parameters for 
the shape and position of the distribution of excited quanta. The numerical results confirm all the features 
obtained from the analytic treatment. 

Under the conditions that guarantee inflation, there are two consecutive but distinct inflationary epochs. 
The first stage features a rapid fall-off of the Hubble parameter and is characterized by the quantum aspects 
of the state. During this first stage the large number of quanta in the excited band are redshifted and build 
up an effective homogeneous classical condensate. The amplitude of this condensate is non-perturbatively 
large, of 0{l/X), as a consequence of the non-perturbatively large number of quanta in the band of excited 
modes. 

The second stage is similar to the classical chaotic scenario and it is driven by the dynamics of this 
effective classical condensate, with vanishing expectation value of the scalar field. Under the tsunami slow-roll 
conditions on the quantum state, the total number of e-folds is more than enough to satisfy the constraints 
of inflationary cosmology. The band of excited wave-vectors if not initially outside the causal horizon, 
becomes superhorizon during the first inflationary stage, therefore these excited states do not modify the 
power spectrum of scalar density perturbations on wavelengths that are of cosmological relevance today. 

Therefore, these tsunami-wave quantum states provide a quantum field theoretical justification of chaotic 
(or in general large field) infiationary models and yield to a microscopic understanding of the emergence of 
classical homogeneous field configurations of large amplitude as an effective collective mode built from the 
large number of quanta in the excited band. 

In addition, we recall that it is necessary to choose an initial state that breaks the $ ^ — $ symmetry in 
classical chaotic scenarios ||, |[ D ■ This is not the case here. We have infiation with zero expectation value 
of the scalar field. 

For completeness we have also studied more general states and established the important difference between 
tsunami (pure or mixed) quantum states leading to inflation, and thermal mixed states which do not lead 
to inflation. 
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APPENDIX A: EQUATIONS OF MOTION IN THE LARGE iV LIMIT AND INITIAL STATES. 



In this appendix we obtain the equations of motion in conformal time for the generalized case in which 
the initial state is determined by a density matrix. The evolution of the functional density matrix is given 
by the Liouville equation in conformal time 



.dp 



H 



5* 



4' 



H 



(Al) 



where the Hamiltonian H is given by eq.(2.1[) to leading order in the large N limit. Consistently with the 
fact that in the large limit the Hamiltonian describes a collection of harmonic oscillators, we propose a 
Gaussian density matrix 



=AA(r)J]exp|- 



TTfc • 7r_fc 



(A2) 
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The hermiticity condition = p for the density matrix impose that Bk must be real. In addition, since 
■k('k,T) is a real field, its Fourier components must obey the hermiticity condition 7f_fe(T) = 7r^(T); thus, 
we can assume A_fc(T) = Ak{T) without loss of generality. 

The evolution equations for Ak{T), M{T) and Bk{T) are obtained from the Liouville eq. (Al) where the 
hamiltonian is given by eq. (2.15). We find 



Al-Bl-a'{T)Lol{T) 



iB, = Bk {Ak - Al) 



(A3) 



where the prime denotes derivative with respect to conformal time T. 

The normalization factor for mixed states Afp{T) is related with the normalization factor of pure states 
Mr{T) by 



where 



AfriT) ==7Vt(0) exp <^ ~i 



dT' 



NVha{T') - 



5](^...^_..)(T')) +y^A,(r' 



8A^ 



Writing Ak in terms of its real and imaginary parts Ak = An ^ + iAj^k, we find that Bk/An^k is a conserved 
quantity. Thus, we can introduce without loss of generality the variables Ar k{T)^ Ai k(^) and 9^ defined 

by 



^fl,fc(T) = Ar.AT) cothGfc 

ArAT) 



AjMT) = Ai,k{T) 



Bk{T) 



sinh 6fc 



where Qk is a time independent real function. 
Introducing the complex variable 

Ak = AR,k + iAi^k 
we see that it obeys the following Ricatti equation 

This equation can be linearized defining 



Ak{r) 



(A4) 

(A5) 
(A6) 

(A7) 



Then eq. (A6) implies that the mode functions Lpk obey 



a\r) 



n{T) 



6 



(A8) 



where Tl{T) is the Ricci scalar. 

The relation ( |A7] ) defines the mode functions ipk{T) up to an arbitrary multiplicative constant that we 
choose such that the wronskian takes the value, 



ipk (fk - fk fk ^ 2i 



(A9) 
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For this choice of the Wronskian the definition (A7) becomes 

1 i d 



2 dr 



In \ipk\ 



(AlO) 



The mass term in eq.(A8) given by eq. (2.17) requires the self-consistent expectation value 



N 

(TTfe • 7r_fc)p 



d^k 



(27r)3 



2\A 



R,k 



Bk 



2A 



coth 



R,k 



Ok 
2 



IVJfeP coth 



2 



Thus, the evolution equations in terms of the mode functions are given by eq. ([A8|) with 



d^k 
(2^ 



coth ( ^ 



(All) 



(A12) 



The evolution equation of the mode functions ipk is the same as the Heisenberg equations of motion for the 
fields, hence we can write the Heisenberg field operators as 



7r(x, T) 



d^k 



x/2(27r)3 



-i fc-x 



(A13) 



Thus, the definit ion ( A7) gives the relation between Schrodinger and Heisenberg pictures, since the functional 
density matrix (A2) is in the Schrodinger picture. 

The expectation value iJj{T) [see eq. (2.13)] in conformal time obeys the following equation of motionQ 



{r)+a\T) 



7^(r) 



^(T) = 



m = ^o ; V''(0) = ^i 



(A14) 
(A15) 



Hence, the evolution equations are given by (A8), (A12) and (A14) with (A15). 

The initial density matrix in the Schrodinguer picture is determined by specifying the initial values of 
AR,k, Ai^k and Ofc. We will take a(0) — 1 and parameterize the initial value of Ak as follows, 



AR.kiO) = ^k , AiM{0)^LOk{0)Sk 



(A16) 



The corresponding initial con diti ons for the mode functions are obtained from eq.(A16) using eq. ( [AT] ) and 
the Wronskian constraint eq.([A9|). These are given by 



<^fc(0) 



ifikiO) = -[ujkiO) Sk + iftk] cpkiO) 



(A17) 



Defin ing the number of particles in terms of the adiabatic eigenstates of the Hamiltonian ( ^.15 ) as in eq. 
(2.27), it is straightforward to find that the initial occupation numbers are given by 



nfc(O) = (rifc(0))p(o) = 



4 ujk ^k 



coth ( ^ 1 - - 



(A18) 



For any mixing parameter Qk 7^ the density matrix represents a mixed state since Bk 7^ 0, a pure initial 
state is obtained by taking Ok = c», in which case Bk — > and the density matrix becomes a product of a 
wave functional times its complex conjugate. 

It is convenient to pass to comoving time, this is achieved by the rescaling of the fields 



^iT{t)) = <t>{t)a{t) , ipki'Tit)) ^ fk{t) a{t) 



(A19) 
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in terms of which the equations of motion are 

4>it) +iH{t) (/)(<) + M^it)<l){t) ^ 



h{t)+ZH{t)h{t) 



M\t) 



fk{t) - 



(2,)3 IA-WPcoth(:^ 



2 



(A20) 



where the dots denote derivative with respect to the comoving time t. The initial conditions for the order 
parameter are its initial value 0(0), and its initial derivative For a(0) = 1, the initial conditions for 

the fluctuations are given by 6^ and 



/fc(0) 



1 



/fc(0) - -[t^fe(O) 4 + H{Q) + ink] /fc(0) 



(A21) 



[Those are the transformed of the initial conditions in conformal time eq. (A17).] 
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